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I.  Summary 

Electrical  Surge  Arrestors  (ESA's)  are  widely  used 
to  protect  susceptible  electrical  and  electronic  components 
sub-systems  and  systems  from  the  potentially  damaging 
effects  of  lightning  static  discharge.  Electromagnetic- 
interference  and  the  Electromagnetic  pulse  arising  from 
potential  nuclear  weapon  detonations.  ESA's  are  very  non- 
linear devices  when  operating  near  or  above  the  DC  break- 
down potential.  The  ESA  model  to  be  described  accounts 
for  a significant  number  of  nonlinear  effects  including  (1) 
streamer  formation.  (2)  plasma  conductivity-.  (3)  glow 
region.  <4  i arc  extinguishing.  (5)  high  frequency  oscilla- 
tions. (fi l Thermionic  effects,  and  (7)  heating  and  heat 
dissipation  to  the  surrounding  medium.  The  model  is  illus- 
trated in  Figure  1 with  a correspondingTlpfinltion  for  each 
parameter  in  the  model. 


Definition  of  Model  Parameters  and  Equations 
1.  Linear  Section 

L Lead  Inductance  (usually  in  nano  henries) 


pL 


Flux  loss  ("Q">  associated  with  L (usually  in 
K ohms) 


Cj,  Stray  Capacitance  (reflected  capacitance  at  the 
ESA  terminals  from  all  other  sources,  leads, 
etc.  ) 

Cf  Gap  Capacitance  (measured  or  calculated  for 
the  Gap) 

2.  Nonlinear  Section 

(a)  R A nonlinear  resistor  dependent  on  the 
* energy  being  dissipated  within  the  gap 

and  transferred  to  the  surrounding  medium 


Equations  : 
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Figure  1.  Electrical  Surge  Arrestor  (ESA)  Model 


II.  Modal  Parameters  and  Equations 

Figure  1 illustrates  the  equivalent  circuit  for  the 
Electrical  Surge  Arrestor  model.  The  linear  portion  of  the 

model  Is  defined  by  the  L,  R . . C and  C elements.  The 
’ pL  g s 

nonlinear  portion  of  the  model  is  also  illustrated  and  basi- 
cally represents  the  complex  reaction  of  initial  streamer/ 
arc  formation  followed  by  plasma  formation  and  Includes  the 
Therm  Ionic  potential  observed  during  the  ON  or  conduction 
phage  of  ESA  operation.  The  various  model  parameters, 
nonlinear  equations  which  utilize  the  parameters  will  be 
defined  and  Illustrated  by  model  application  and  comparison 
to  test  results.  The  various  test  circuits  are  also  shown  for 
the  benefit  of  other  Investigators  wishing  to  characterize 
ESA's  or  spark  gaps  in  a manner  described  herein. 


whe  re : 
R 


(whichever  Is  less) 

Instantam-ous  Resistance  of  the  gap  lohm) 
Instantaneous  Voltage  across  the  gap  (voltst 


1^  Instantaneous  Current  through  the  gap  (amps) 

W Instantaneous  Energy  in  the  gap  (Joules) 

W)N  Instantaneous  Energy  generated  as  Input  to 
the  gap  (Joules) 

W„t,_  * Instantaneous  Energy  removed  from  the  gap 
to  the  surrounding  medium  (joules! 


Scale  factor  relating  R to  W (tl-Joules) 


» 

r 


i 


I 


i 


Pmj  Maximum  Instantaneous  power  lost  to  the 
surrounding  medium  (watts) 

r . Dissipation  time  constant  urn  average 
1 thermal  time  constant  for  establishing  the 
rate  at  which  heat  energy  (joules)  is  dissi- 
pated to  the  surrounding  medium)  (sec) 

R Maximum  resistance  of  the  air  gap  (ohms) 

gmax 

Rgmin  Minimum  resistance  of  the  air  gap  (ohms) 
"min  Minimum  energy  in  the  air  gap  (joules) 

V Gap  breakdown  threshold  voltage  (volts) 
UU  V'DB  = lOOKV/inch  in  air  at  STP 


t Time  for  Streamer  Formation  (r 

or  or 

increases  with  gap  width) 


I 


gmin 


Instantaneous  value  of  arc  current  below 
which  the  arc  extinguishes. 


The  following  controls  are  employed: 


No.  1 If  V » and  W » 0.  then  initiate  timer 

at  t (time  at  which  the  condition  for  streamer 

o 

formation  has  been  achieved) 

No.  2 At  t tQ  4 t initiate  computation  of  W 

and  subsequent  modification  of  R (see 
Figure  2 characteristic  curve)  B 

No.  3 Continue  to  monitor  W until  W s W . . Also 

min 

monitor  R and  lim‘t  R to  R bv  com- 
g g gmin 

paring  R <W)  to  R 
v * g gmin 


P 

< 

x 


S 


STREAM! H 


V (OOV'OIV 


(1)1  LOW  fRtQUtNCY  I VSV  CMARACTi  HISTICS  Of  lal 


No.  4 Return  R to  R when  gap  is  fully 
g gmax 

extinguished  I < I 

g gmin 


Figure  2.  Non-Linear  Section  of  FSA  Model  and 
Low  Frequency  I (V)  Characteristics 


(b)  Thermionic  Kffect  The  phenomenon  that 
accounts  for  the  high  ON  voltages  observed  in 
ESA’s  (over  100  volts),  whereby  the  emitting  ESA 
block  effectively  becomes  a cold  cathode,  while  the 
collecting  ESA  block  becomes  the  plate,  in  what 
Is  essentially  Thermionic  reaction,  namely,  the 
plasma  creates  a junction  which  is  represented 
mathematically  by  the  following  equations: 


RFl 

lgl  (exp  VR/MjS)  - 1) 

(5) 

RF2 

(exp  -VR/M2»)  - 1) 

(6) 

where: 

V„  Thermionic  junction  potential  (volts) 

K 

l_  , * A pseudo- saturation  current  for  the 
‘ Thermionic  Rectifier  (and  are  also  func- 

tions of  trmperature  and  effective  plasma 
area) 

M,,  M„  Multiplier  for  empirical  fit  (M  has  a range 
of  50  to  500  dopenillng  on  the  particular 
ESA)  (non-dimensional) 


In  addition,  there  Is  a slight  tendency  for  these 
"Thermionic  Rectifiers"  to  "store"  charge  similar  to  a 
semiconductor  rectifier,  hence,  the  net  current  in  each 
Thermionic  Rectifier  is  formulated  as  follows: 


HI 

'hit 

♦ r I 

It  RFl 

(71 

I ~ 

4 t i 

(8) 

H2 

RF2 

K RF2 

where: 

t Thermionic  electron  recombination  time  (which 
is  believed  to  be  on  the  order  of  a few 
nanoseconds ) 

0 kT/q  = . 026  (a  27°C  (T  is  set  & 300°K)  (9) 

A set  of  typical  ESA  model  data  is  shown  in  Table  1. 

Using  non-linear  model  shown  in  Figure  Land  the  data 
shown  in  Table  I.  it  was  possible  to  analytically  obtain 
1 versus  V characteristics  showing  the  extremely  non-linear 
behavior  of  the  gap  for  (1)  Streamer,  (2)  Arc/plasma,  and 
(3)  Thermionic  formation  including  extinguishing  of  the  arc. 
This  is  shown  in  Figure  2. 
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Table  l.  F.xample  of  ESA  Data 


V - 2(|w) 


Parameter 

Value 

Units 

L 

130 

nh 

Linear 

V 

2500 

ohms 

Section 

Cg 

3.  18 

pfd 

C 

S 

1 

pfd 

Ki 

2/'3 

ohm -Joules 

Pml 

1500 

Watts 

Td 

.1 

usee 

R 

gmax 

1/+12 

ohms 

Nonlinear  „ 
Section 

R 

gmtn 

. 1 

ohms 

W . 
nun 

/-io 

Joules 

VDB 

185 

Y'olts 

tSF 

1 

nsec 

'gmtn 

5 

ma 

1 

8 

1 

ua 

Thermionic 

M 

500 

- 

H 

. 1 

nsec 

III  ESA  CHARACTERIZATION 
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(a)  SMALL  SIGNAL  EQUIVALENT  CIRCUIT 
(FIGURE  1 LINEAR  SECTIONI 


HOR  ■ 1 MHz  THROUGH  100  MHz 

(1)1  LABORATORY  TEST  RESULTS 
(USING  HP  NETWORK  ANALYSER) 


Several  types  of  KSA’s  were  modeled  and 
characterized  (Audio,  Power  and  Antenna  h:S/Vsi  for 
both  the  linear  and  non-linear  components  of  the 
model. 


Figure  3.  Small  Signal  Equivalent  Circuit  for 
Audio  t'SA  (Below  Breakdown ) and 
laboratory  Test  Response 


"Hie  linear  characteristics  can  be  obtained  quite  readily 
bv  using  a network  analyzer  and  solving  for  the  values  of  L, 

c'  , C , and  R , . The  result  for  an  Audio  ESA  is  shown  in 
S g pL 

Figure  3, 

The  nonlinear  characteristics  are  considerably  more 
difficult  to  obtain. 

The  following  tests  were  conducted  to  obtain  the  non- 
linear characteristics: 

1.  Sawtooth  Oscillator  Test 

The  first  test  that  can  be  performed  is  Rhown  in  Fig- 
ure 4 where  the  ESA  is  used  as  the  nonlinear  element  which 
produces  a sawtooth  oscillation.  Resistor.  Rj.  must  be 

sufficiently  large  so  that  the  'holding  current'  wi!'  not  sus- 
tain a very  low  electron  leakage  and  prevent  oscillation. 

This  determines  one  of  the  critical  parameters  namely, 
the  minimum  power  to  sustain  the  arc.  This  may  not  he  the 
minimum  power  required  to  completely  describe  the  inter- 
action of  the  plasma  with  the  surrounding  medium  however. 
The  dynamic  resistance  of  the  arc  can  appear  quite  high 
even  though  the  capacitor.  Cj  is  a very  low  impedance 

source  at  the  switch  point  and  hence,  would  be  expected  to 
discharge  rapidly  giving  rise  to  a high  current  pulse  through 
the  ESA;  this  is  not  the  case  for  this  or  other  ESA's,  the 
reason  for  this  apparent  high  discharge  impedance  in  this 
sawtooth  oscillator  configuration  Is  postulated  to  be  due  to 


the  relatively  high  resistance  of  R and  the  relatively  high 
Imiwdance  of  the  onset  of  the  Thermionic  conduction 


'gap 


2.  Fast  Rise  Time  Test 


When  the  ESA  is  subject  to  a high  voltage,  fast  rise- 
time  input,  two  characteristics  become  apparent.  E'irst. 
the  firing  |K>tential  increases  and  second,  the  ON  voltage  is 
only  somewhat  higher  than  during  the  sawtooth  operation. 

The  dynamic  resistance  becomes  much  lower  which  satisfies 

the  functional  relationship  between  R and  Mo/I  (both 

gap  B gap 

resistors  being  inverse  to  arc  current  density).  The  high 
dv/dt  input  results  in  a voltage  breakdown  vs  rise-time 
characteristics  as  shown  in  Figure  5.  There  are  two 
regions  of  the  curve  which  are  worthy  of  discussion. 

Region  I shows  a fairly  gradual  increase  in  apparent  gap 
breakdown  with  Increasing  dv/dt.  This  increase  is  attribut- 
able to  the  Interaction  of  the  arc  formation  and  energy  dissi- 
pation coupled  with  the  effect  of  the  reactive  linear  elements 
(L's  and  C's)  of  the  ESA.  in  Region  II  the  apparent  break- 
down of  the  ESA  vs  dv/dt  Increased  more  rapidly.  This  is 
attributed  to  the  time  required  for  streamer  formation. 

This  occurs  prior  to  arc  formation.  The  net  effect  of  this 
ESA  response  to  very  high  values  of  dv/dt  (or  high  fre- 
quency) reduces  ESA  effectiveness  for  the  high  frequency 
components  of  disturbance. 
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(a)  FAST  RISE  TIME  TEST  CIRCUIT  (R,  1 MEG.C,  • 0 02UFD. 

R7.R3. 5I0W.  C,  VARIED  FROM  50  PFD  THROUGH  350  PFD) 
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Figure  4.  Sawtooth  Oaclllator  Teat  ResultB  and 
Comparlaon  to  Computer  Model  Response 
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> REGION  I - 


0.01  0.1  1.0  10  100 
RISE  TIME  Idv/dtl  (KV/nSECI 

Id)  PLOT  OF  INCREASED  BREAKDOWN  VOLTAGE  V0  (APPARENT) 
VS  RISE  TIME  (dv/dt)  OR  EQUIVALENT  FREQUENCY  I 


Figure  5.  Characterization  of  Increase  in  Apparent 
Voltage  Breakdown  (or  Overshoot)  vs  Rise  Time 
Test  Circuit.  Results  and  Plot  of 
Vg  vs  dv/dt 
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3.  Damped  Sint'  Wave  Test 

In  practice  a damped  sinewave  is  frequently  used  to 
characterize  component  and/or  subsystem  response  to 
Klectrical  Surges.  As  a result,  the  ESA  and  corresponding 
model  were  also  characterized  using  various  damped  sine- 
wave  input  stimulus.  The  result  of  this  test  and  model 
simulation  as  illustrated  in  Figure  6.  Note  that  the  ESA's 
can  open  after  the  initial  shorting  if  the  sustaining  energy 
is  insufficient  to  cause  the  arc  to  sustain.  This  is  particu- 
larly true  for  lower  frequencies  (~1  Mllz  and  belowi  and 
higher  impedance  n K'.l)  sources. 

IV.  Example  Computer  Codmg  for  the  ESA 
Models  Using  SYSCAP  II 

The  SYSCAP  II  program  was  used  to  perform  several 
studies  involving  the  spark  gaps  and  surge  arrestors  dis- 
cussed in  Section  III.  This  section  will  illustrate  how  some 
of  these  models  were  programmed  using  the  SYSCAP  II 
ALCAP  and  TRACAP  subprograms. 

ALCAP:  The  ALCAP  (A.  C.  Linear  Computer  Analy- 
sis Program!  was  used  to  model  the  linear  network  shown 
in  Figure  3.  Figure  7 illustrates  the  coding  list  used 
to  obtain  the  results  shown  in  Figure  3.  In  addition  to  the 
amplitude  and  phase  vs  frequency  response,  the  ALCAP 
program  was  also  used  to  determine  the  sensitivity  of  the 
amplitude  and  phase  response  to  component  tolerance  at 
various  frequencies.  This  capability  of  the  program  was 
very  useful  in  characterizing  the  networks  since  it  identified 
which  elements  were  the  most  critical.  It  showed  that  L 
and  C were  very  critical  around  resonance  while  R . and 

g PL 

C were  not  as  sensitive. 

8 

TRACAP:  The  TRACAP  (TRAnsient  Computer  Anal- 
ysis Program)  was  used  to  model  the  combined  linear  and 
non-lir.ear  model  shown  in  Figure  I.  Figure  8 illustrates 
an  example  of  TRACAP  coding  used  to  obtain  the  non- 
linear !(v)  characteristics  (as  shown  in  Figure  2).  It  should 
be  noted  that  there  Is  a significant  negative  resistance 
region  which  is  characteristic  of  spark  gaps  which  has  been 
very  difficult  to  model  in  the  past  but  posed  no  significant 
problems  using  the  TRACAP  code  due  to  the  conservative 
convergence  algorithms  that  have  been  used. 

Another  example  of  TRACAP  application  for  ESA 
modeling  Is  illustrated  in  Figure  9 where  a damped  sine- 
wave  stimulus  simulates  a potential  surge  disturbance. 

The  results  are  compared  to  test  data  shown  in  Figure  10. 

In  this  case  two  functions  are  used,  namely  FUNC  1 to 
characterize  the  ESA  resistance  IWW)  and  FUNC  2 to 

generate  the  exponentially  damped  sine  wave- 

fg  “ A«  sin  (wt). 

These  examples  have  been  given  so  that  scientists  and 
engineers  who  are  interested  in  this  type  of  model  (surge 
protection  analysis/design)  can  readily  avail  themselves  of 
this  capability  via  the  SYSCAP  program. 


V.  Conclusions 

The  ESA  model  presented  in  this  paper  is  relatively 
easy  to  use  (In  an  appropriate  computer  program  such  as 
SYSCAP  II)  and  permits  a detailed  evaluation  spark  gap 
protection  to  a variety  of  electrical  surge  stimulus.  The 
model  accounts  for  the  extremely  nonlinear  behavior  of  the 
ESA  gap  Including  intermittent  firing,  variation  in  ON 
Impedance  and  other  related  phenomenon.  Examples  of 
model  application  demonstrate  the  behavior  of  the  model 
in  the  linear  and  non-llnear  regions. 


EXAMPLE  OF  PULSES  MODEL  IC/RI  USING  SIMPLE 
EQUIVALENT  CIRCUIT  LAB  TEST  RESULTS 
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0 5pS/DIV 
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It)  ESA  MODEL  RESPONSE  (INCLUDING  SHOR T 

INDUCTANCE  LOOP  RINGING)  (TEKTRONIX  *0)0) 


Figure  6.  Comparison  of  ESA  lab  Tests 
(Damped  Sinewave)  and  ESA  Model 
Response  iTektronix  4010  Display) 
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Figure  7.  SYSCAP  II  (ALCAP)  Coding  and 
Amplltude/Phase  Response  Using  a TI 
Silent  700  Terminal 


VI. 


SYSCAP  II  User  Information  Manual  (Publication 
No.  76070600).  Available  through  CDC  Cybernet 
Data  Services  Publications.  P. O.  Box  O.  HQW05F, 
Minneapolis.  MN  55440 
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Figure  8.  Illustration  of  Interactive  SYSCAP  II  (TRACAP) 
Coding  and  Response  Using  a TI  700  Terminal 
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Figure  9.  Illustration  of  Interactive  SY8CAP  II 
(TRACAP)  Coding  and  Response  for  a 
Damped  SInewave  Input 
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Figure  10.  Test  Circuit  and  Response  Comparison 
for  a Damped  SInewave 
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